Molecular characterization of the cytochrome P450 enzyme CYP18A1 in Henosepilachna vigintioctopunctata.

In insects, the expression of 20E response genes that initiate metamorphosis is triggered by a pulse of 20-hydroxyecdysone (20E). The 20E pulse is generated through two processes: synthesis, which increases its level, and inactivation, which decreases its titer. CYP18A1 functions as an ecdysteroid 26-hydroxylase and plays a role in 20E removal in several representative insects. However, applying 20E degradation activity of CYP18A1 to other insects remains a significant challenge. In this study, we discovered high levels of Hvcyp18a1 during the larval and late pupal stages, particularly in the larval epidermis and fat body of Henosepilachna vigintioctopunctata, a damaging Coleopteran pest of potatoes. RNA interference (RNAi) targeting Hvcyp18a1 disrupted the pupation. Approximately 75% of the Hvcyp18a1 RNAi larvae experienced developmental arrest and remained as stunted prepupae. Subsequently, they gradually turned black and eventually died. Among the Hvcyp18a1-depleted animals that successfully pupated, around half became malformed pupae with swollen elytra and hindwings. The emerged adults from these deformed pupae appeared misshapen, with shriveled elytra and hindwings, and were wrapped in the pupal exuviae. Furthermore, RNAi of Hvcyp18a1 increased the expression of a 20E receptor gene (HvEcR) and four 20E response transcripts (HvE75, HvHR3, HvBrC, and HvαFTZ-F1), while decreased the transcription of HvßFTZ-F1. Our findings confirm the vital role of CYP18A1 in the pupation, potentially involved in the degradation of 20E in H. vigintioctopunctata.

Six complete mitochondrial genomes of ground beetles from the Harpalinae and Carabinae (Coleoptera, Carabidae) with phylogenetic analysis based on mitogenomic data.

In this study, we employed high-throughput sequencing technology to determine the complete mitochondrial genomes of six ground beetles, encompassing five Harpalinae species and one Carabinae species. The sizes of mitochondrial genomes ranged from 15,334 to 16,972 bp, encompassing 37 genes, including 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. Furthermore, each species was found to possess a putative control region. Combining with 65 published mitochondrial genome sequences of Carabidae as ingroups and four species from Trachypachidae, Gyrinidae and Dytiscidae as outgroups, we conducted phylogenetic analyses utilizing Maximum likelihood and Bayesian inference methods. Moreover, we reconstructed a species tree of Carabidae based on mitochondrial genome data using the coalescent-based species tree method (ASTRAL). The results revealed that the family Carabidae was not a monophyletic group. The subfamily Harpalinae was supported to be a monophyletic group in Maximum likelihood analysis. Although the subfamily Carabinae was found to be nonmonophyletic in the concatenation analyses under both Maximum likelihood and Bayesian inference criteria, it was identified as a monophyletic group in the species tree analysis.

The first chromosome-level genome of the stag beetle Dorcus hopei Saunders, 1854 (Coleoptera: Lucanidae).

Stag beetles (Coleoptera: Lucanidae) represent a significant saproxylic assemblage in forest ecosystems and are noted for their enlarged mandibles and male polymorphism. Despite their relevance as ideal models for the study of exaggerated mandibles that aid in attracting mates, the regulatory mechanisms associated with these traits remain understudied, and restricted by the lack of high-quality reference genomes for stag beetles. To address this limitation, we successfully assembled the first chromosome-level genome of a representative species Dorcus hopei. The genome was 496.58 Mb in length, with a scaffold N50 size of 54.61 Mb, BUSCO values of 99.8%, and 96.8% of scaffolds anchored to nine pairs of chromosomes. We identified 285.27 Mb (57.45%) of repeat sequences and annotated 11,231 protein-coding genes. This genome will be a valuable resource for further understanding the evolution and ecology of stag beetles, and provides a basis for studying the mechanisms of exaggerated mandibles through comparative analysis.

Vertically inherited microbiota and environment modifying behaviours conceal genetic variation in dung beetle life history.

Diverse organisms actively manipulate their (sym)biotic and physical environment in ways that feed back on their own development. However, the degree to which these processes affect microevolution remains poorly understood. The gazelle dung beetle both physically modifies its ontogenetic environment and structures its biotic interactions through vertical symbiont transmission. By experimentally eliminating (i) physical environmental modifications and (ii) the vertical inheritance of microbes, we assess how environment modifying behaviour and microbiome transmission shape heritable variation and evolutionary potential. We found that depriving larvae of symbionts and environment modifying behaviours increased additive genetic variance and heritability for development time but not body size. This suggests that larvae's ability to manipulate their environment has the potential to modify heritable variation and to facilitate the accumulation of cryptic genetic variation. This cryptic variation may become released and selectable when organisms encounter environments that are less amenable to organismal manipulation or restructuring. Our findings also suggest that intact microbiomes, which are commonly thought to increase genetic variation of their hosts, may instead reduce and conceal heritable variation. More broadly, our findings highlight that the ability of organisms to actively manipulate their environment may affect the potential of populations to evolve when encountering novel, stressful conditions.

Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles.

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

Comparative transcriptomic and metabolomics analysis of modified atmosphere responses in Tribolium castaneum (Coleoptera: Tenebrionidae).

Modified atmosphere is effective in controlling Tribolium castaneum Herbst, but it has adaptations. Comprehending the potential mechanism of resistance to T. castaneum in a modified atmosphere will help advance related management methods. This study conducted a comparative transcriptomic and metabolomic analysis to understand the physiological mechanism of T. castaneum in adapting to CO2 stress. Results showed that there were a large number of differentially expressed genes (DEGs) in T. castaneum treated with different concentrations of CO2. Gene ontology (GO) analysis revealed significant enrichment of DEGs mainly in binding, catalytic activity, cell, membrane, membrane part, protein-containing complex, biological regulation, and cellular and metabolic process. Kyoto Encyclopedia of Genes and Genomes analysis showed that different treatments had different effects on the metabolic pathways of T. castaneum. DEGs induced by 25% CO2 were involved in arginine and proline metabolism, and 50% air + 50% CO2 treatment affected most kinds of metabolic pathways, mainly the signal transduction pathway, including PI3K-Akt signaling pathway, AMPK signaling pathway, neurotrophin signaling pathway, insulin signaling pathway, and thyroid hormone signaling. Ribosome and DNA replication were enriched under high CO2 stress (75% and 95%). The metabolomics revealed that different concentrations of CO2 treatments might inhibit the growth of T. castaneum through acidosis, or they may adapt to anoxic conditions through histamine and N-acetylhistamine. Multiple analyses have shown significant changes in histamine and N-acetylhistamine levels, as well as their associated genes, with increasing CO2 concentration. In conclusion, this study comprehensively revealed the molecular mechanism of T. castaneum responding to CO2 stress and provided the basis for an effectively modified atmosphere in the T. castaneum.

The construction and evaluation of secretory expression engineering bacteria for the trans-Cry3Aa-T-HasA fusion protein against the Monochamus alternatus vector.

Pine wood nematode disease is currently the most deadly forest disease in China, and the Monochamus alternatus is its primary vector. Controlling the M. alternatus is crucial for managing pine wood nematode disease. This study, based on the selected HasA (pGHKW4) secretory expression vector, used electroporation to combine the genetically modified high-toxicity toxin Cry3Aa-T with the entomopathogenic bacterium Yersinia entomophaga isolated from the gut of the M. alternatus. The SDS-PAGE and Western blotting techniques were employed to confirm the toxin protein's secretion capability. The engineered bacteria's genetic stability and effectiveness in controlling M. alternatus were assessed for their insecticidal activity. The results of the SDS-PAGE and Western blotting analyses indicate that the HasA system effectively expresses toxin protein secretion, demonstrates certain genetic stability, and exhibits high insecticidal activity against M. alternatus. This study constructed a highly toxic entomopathogenic engineered bacterial strain against M. alternatus larvae, which holds significant implications for controlling M. alternatus, laying the foundation for subsequent research and application of this strain.

Resistance management and integrated pest management insights from deployment of a Cry3Bb1+ Gpp34Ab1/Tpp35Ab1 pyramid in a resistant western corn rootworm landscape.

In Nebraska USA, many populations of western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte, now exhibit some level of resistance to all corn rootworm-active Bacillus thuringiensis Berliner (Bt) proteins expressed in commercial hybrids. Therefore, a study was conducted in northeast Nebraska from 2020-2022 to reevaluate current corn rootworm management options in continuous maize (consecutive planting for ≥2 years). Results from on-farm experiments to evaluate a standard soil-applied insecticide (Aztec® 4.67G) in combination with non-rootworm Bt or rootworm-active Bt pyramided maize (Cry3Bb1 + Gpp34Ab1/Tpp35Ab1) are reported within the context of WCR Bt resistance levels present. Corrected survival from Bt pyramid single-plant bioassays (<0.3, 0.3-0.49, >0.5) was used to place populations into 3 resistance categories. Variables evaluated included root injury, adult emergence, proportion lodged maize, and grain yield. Key results: A composite analysis of all populations across resistance levels indicated that addition of soil insecticide to Bt pyramid significantly reduced adult emergence and lodging but did not significantly increase root protection or yield. Within and among resistance category analyses of root injury revealed that the Bt pyramid remained highly efficacious at any non-rootworm Bt root injury level when resistance was absent or low. When corrected survival was >0.3, mean Bt pyramid root injury tracked more closely in a positive linear fashion with mean non-rootworm Bt root injury (rootworm density x level of resistance interaction). Similar trends were obtained for adult emergence but not yield. Mean Bt pyramid root injury rating was <0.75 in most populations with Bt resistance, which contributed to no significant yield differences among categories. Results are discussed within the context of IPM:IRM tradeoffs and the need to reduce WCR densities in this system to decrease the impact of the density x resistance interaction to bridge use of current pyramids with new technologies introduced over the next decade.

Three chemosensory proteins enriched in antennae and tarsi of Rhaphuma horsfieldi differentially contribute to the binding of insecticides.

Antennae and legs (primarily the tarsal segments) of insects are the foremost sensory organs that contact a diverse range of toxic chemicals including insecticides. Binding proteins expressed in the two tissues are potential molecular candidates serving as the binding and sequestering of insecticides, like chemosensory proteins (CSPs). Insect CSPs endowed with multiple roles have been suggested to participate in insecticide resistance, focusing mainly on moths, aphids and mosquitos. Yet, the molecular underpinnings underlying the interactions of cerambycid CSPs and insecticides remain unexplored. Here, we present binding properties of three antenna- and tarsus-enriched RhorCSPs (RhorCSP1, CSP2 and CSP3) in Rhaphuma horsfieldi to eight insecticide classes totaling 15 chemicals. From the transcriptome of this beetle, totally 16 CSP-coding genes were found, with seven full-length sequences. In phylogeny, these RhorCSPs were distributed dispersedly in different clades. Expression profiles revealed the abundant expression of RhorCSP1, CSP2 and CSP3 in antennae and tarsi, thus as representatives for studying the protein-insecticide interactions. Binding assays showed that the three RhorCSPs were tuned differentially to insecticides but exhibited the highest affinities with hexaflumuron, chlorpyrifos and rotenone (dissociation constants <13 µM). In particular, RhorCSP3 could interact strongly with 10 of tested insecticides, of which four residues (Tyr25, Phe42, Val65 and Phe68) contributed significantly to the binding of six, four, three and four ligands, respectively. Of these, the binding of four mutated RhorCSP3s to a botanical insecticide rotenone was significantly weakened compared to the wildtype protein. Furthermore, we also evidenced that RhorCSP3 was a broadly-tuned carrier protein in response to a wide variety of plant odorants outside insecticides. Altogether, our findings shed light on different binding mechanisms and odorant-tuning profiles of three RhorCSPs in R. horsfieldi and identify key residues of the RhorCSP3-insecticide interactions.

Cloning and characterization of luciferase from an Asian firefly Pygoluciola qingyu and its comparison with other beetle luciferases.

The bioluminescence system of luminescent beetles has extensive applications in biological imaging, protein labeling and drug screening. To explore wild luciferases with excellent catalytic activity and thermal stability, we cloned the luciferase of Pygoluciola qingyu, one species living in areas of high temperature and with strong bioluminescence, by combining transcriptomic sequencing and reverse transcription polymerase chain reaction (RT-PCR). The total length of luciferase gene is 1638 bp and the luciferase consists 544 amino acids. The recombinant P. qingyu luciferase was produced in vitro and its characteristics were compared with those of eight luciferases from China firefly species and two commercial luciferases. Compared with these luciferases, the P. qingyu luciferase shows the highest luminescence activity at room temperature (about 25-28 â„ƒ) with similar KM value for D-luciferin and ATP to the Photinus pyralis luciferase. The P. qingyu luciferase activity was highest at 35 â„ƒ and can keep high activity at 30-40 â„ƒ, which suggests the potential of P. qingyu luciferase for in vivo and cell application. Our results provide new insights into P. qingyu luciferase and give a new resource for the application of luciferases.

Coconut rhinoceros beetle digestive symbiosis with potential plant cell wall degrading microbes.

Coconut rhinoceros beetle (CRB, Oryctes rhinoceros) is an invasive palm pest whose larvae eat wood, yet lack the necessary digestive enzymes. This study confirmed endogenous CRB cellulase is inactive, suggesting microbial fermentation. The inner lining of the CRB hindgut has tree-like structures covered with a conspicuous biofilm. To identify possible symbionts, 16 S rRNA amplicon sequencing was used on individuals from across Taiwan. Several taxa of Clostridia, an anaerobic class including many cellulolytic bacteria, were highly abundant in most individuals from all locations. Whole metagenome sequencing further confirmed many lignocellulose degrading enzymes are derived from these taxa. Analyses of eggs, larvae, adults, and soil found these cellulolytic microbes are not transmitted vertically or transstadially. The core microbiomes of the larval CRB are likely acquired and enriched from the environment with each molt, and enable efficient digestion of wood.

Mammal dung-dung beetle trophic networks: an improved method based on gut-content DNA.

Background: Dung beetles provide many important ecosystem services, including dung decomposition, pathogen control, soil aeration, and secondary seed dispersal. Yet, the biology of most dung beetles remains unknown. Natural diets are poorly studied, partly because previous research has focused on choice or attraction experiments using few, easily accessible dung types from zoo animals, farm animals, or humans. This way, many links within natural food webs have certainly been missed. In this work, we aimed to establish a protocol to analyze the natural diets of dung beetles using DNA gut barcoding. Methods: First, the feasibility of gut-content DNA extraction and amplification of 12s rDNA from six different mammal dung types was tested in the laboratory. We then applied the method to beetles caught in pitfall traps in Ecuador and Germany by using 12s rDNA primers. For a subset of the dung beetles caught in the Ecuador sampling, we also used 16s rDNA primers to see if these would improve the number of species we could identify. We predicted the likelihood of amplifying DNA using gut fullness, DNA concentration, PCR primer, collection method, and beetle species as predictor variables in a dominance analysis. Based on the gut barcodes, we generated a dung beetle-mammal network for both field sites (Ecuador and Germany) and analyzed the levels of network specificity. Results: We successfully amplified mammal DNA from dung beetle gut contents for 128 specimens, which included such prominent species as Panthera onca (jaguar) and Puma concolor (puma). The overall success rate of DNA amplification was 53%. The best predictors for amplification success were gut fullness and DNA concentration, suggesting the success rate can be increased by focusing on beetles with a full gut. The mammal dung-dung beetle networks differed from purely random network models and showed a moderate degree of network specialization (H2': Ecuador = 0.49; Germany = 0.41). Conclusion: We here present a reliable method of extracting and amplifying gut-content DNA from dung beetles. Identifying mammal dung via DNA reference libraries, we created mammal dung-dung beetle trophic networks. This has benefits over previous methods because we inventoried the natural mammal dung resources of dung beetles instead of using artificial mammal baits. Our results revealed higher levels of specialization than expected and more rodent DNA than expected in Germany, suggesting that the presented method provides more detailed insights into mammal dung-dung beetle networks. In addition, the method could have applications for mammal monitoring in many ecosystems.

scRNA-seq Reveals Novel Genetic Pathways and Sex Chromosome Regulation in Tribolium Spermatogenesis.

Spermatogenesis is critical to sexual reproduction yet evolves rapidly in many organisms. High-throughput single-cell transcriptomics promises unparalleled insight into this important process but understanding can be impeded in nonmodel systems by a lack of known genes that can reliably demarcate biologically meaningful cell populations. Tribolium castaneum, the red flour beetle, lacks known markers for spermatogenesis found in insect species like Drosophila melanogaster. Using single-cell sequencing data collected from adult beetle testes, we implement a strategy for elucidating biologically meaningful cell populations by using transient expression stage identification markers, weighted principal component clustering, and SNP-based haploid/diploid phasing. We identify populations that correspond to observable points in sperm differentiation and find species specific markers for each stage. Our results indicate that molecular pathways underlying spermatogenesis in Coleoptera are substantially diverged from those in Diptera. We also show that most genes on the X chromosome experience meiotic sex chromosome inactivation. Temporal expression of Drosophila MSL complex homologs coupled with spatial analysis of potential chromatin entry sites further suggests that the dosage compensation machinery may mediate escape from meiotic sex chromosome inactivation and postmeiotic reactivation of the X chromosome.

De novo assembly and annotation of Popillia japonica's genome with initial clues to its potential as an invasive pest.

BACKGROUND: The spread of Popillia japonica in non-native areas (USA, Canada, the Azores islands, Italy and Switzerland) poses a significant threat to agriculture and horticulture, as well as to endemic floral biodiversity, entailing that appropriate control measures must be taken to reduce its density and limit its further spread. In this context, the availability of a high quality genomic sequence for the species is liable to foster basic research on the ecology and evolution of the species, as well as on possible biotechnologically-oriented and genetically-informed control measures. RESULTS: The genomic sequence presented and described here is an improvement with respect to the available draft sequence in terms of completeness and contiguity, and includes structural and functional annotations. A comparative analysis of gene families of interest, related to the species ecology and potential for polyphagy and adaptability, revealed a contraction of gustatory receptor genes and a paralogous expansion of some subgroups/subfamilies of odorant receptors, ionotropic receptors and cytochrome P450s. CONCLUSIONS: The new genomic sequence as well as the comparative analyses data may provide a clue to explain the staggering invasive potential of the species and may serve to identify targets for potential biotechnological applications aimed at its control.

Phylogenetic analyses of the subfamily Scarabaeinae (Coleoptera: Scarabaeidae) provide new insights into the Mexican Transition Zone theory.

Being areas of biotic overlap located between biogeographic regions, transition zones function as natural laboratories. The present study explores the phylogenetic history of the dung beetle subfamily Scarabaeinae, in order to present an evolutionary scenario that allows inference of the biogeographic history of the Mexican Transition Zone (MTZ) and integration of the distributional patterns of its biota. The species sampling included 94 New World taxa (93 species of Scarabaeinae and one species of Aphodiinae). The phylogenetic relationships of the main clades recovered in our study were supported with PP values 0.95. Based on the BAYAREALIKE model to reconstruct the ancestral distributional patterns of Scarabaeinae, we inferred a complex scenario with 19 dispersal events, 15 vicariance events, and three extinctions. We suggest that the Ancient Neotropical and Tropical Paleoamerican patterns represent the most likely ancestral distributional patterns for the Scarabaeinae of the MTZ, which probably settle there during the Eocene-Oligocene. The rest of the Scarabaeinae distributional patterns were assembled in subsequent periods. The results suggest that the MTZ had two separate formation stages: a Paleo-MTZ (Eocene-Miocene) and a current MTZ (Pliocene-Anthropocene). We conclude that the evolutionary history as well as the dispersal-vicariance scenario for the Scarabaeinae of the MTZ fits the out of the tropics model.

Odorant binding protein 18 increases the pathogen resistance of the imported willow leaf beetle, Plagiodera versicolora.

Background: Insect odorant-binding proteins (OBPs) are a class of small molecular weight soluble proteins. In the past few years, OBPs had been found to work as carriers of ligands and play a crucial role in olfaction and various other physiological processes, like immunity. A subset of insect OBPs had been found to be expressed differently and play a function in immunity of fungal infection. However, there are few studies on the role of OBPs in immunity of bacterial infection. Methods: To identify the immune-related OBPs of Plagiodera versicolora after infected by Pseudomonas aeruginosa, we determined the mortality of P. versicolora to P. aeruginosa and selected the time point of 50% mortality of larvae to collect samples for RNA-seq. RNAi technology was used to investigate the function of immune-related OBPs after P. aeruginosa infection. Results: RNA-seq data shows that PverOBP18 gene significantly up-regulated by 1.8-fold and further RT-qPCR affirmed its expression. Developmental expression profile showed that the expression of PverOBP18 was highest in the pupae, followed by the female adults, and lower in the 1st-3rd larvae and male adults with lowest in eggs. Tissue expression profiling showed that PverOBP18 was dominantly expressed in the epidermis. RNAi knockdown of PverOBP18 significantly reduced the expression of bacterial recognition receptor gene PGRP and antibacterial peptide gene Attacin and reduced the resistance of P. versicolora to P. aeruginosa infection. Conclusion: Our results indicated that PverOBP18 gene increased the pathogen resistance of P. versicolora by cooperating with the immune genes and provided valuable insights into using OBPs as targets to design novel strategies for management of P. versicolora.

Binding Properties of Odorant Binding Protein 37 in Plagiodera versicolora to Host Volatile, o-Cymene.

The chemosensory system plays an important role in the host plants location. Plagiodera versicolora (Coleoptera: Chrysomelidae) is a worldwide leaf-eating forest pest that feeds exclusively on salicaceous trees. There is no function study of odorant binding proteins (OBPs) in P. versicolora. In the current study, we found that PverOBP37 has a high expression in male and female antennae, heads, and legs by quantitative real-time PCR. The binding properties of PverOBP37 to 18 host plant volatiles were determined by fluorescence competition binding assays. The results showed that PverOBP37 could bind to the host plant volatile, o-cymene. Furthermore, four candidate key amino acid residues (F8, Y50, F103, and R107) of PverOBP37 to o-cymene were identified by molecular docking. The functional assay to confirm Y50, F103, and R107 mutations were key amino acid residues of PverOBP37 involved in the binding to o-cymene. Knockdown of PverOBP37 and Y-tube behavioral bioassays of mated females led to a significantly reduced attraction to o-cymene. This study not only revealed the molecular mechanism of PverOBP37 but also suggested that PverOBP37 is essential to detect host plant volatiles as cues to search for egg-laying sites in P. versicolora.

Species delimitation, discovery and conservation in a tiger beetle species complex despite discordant genetic data.

In an age of species declines, delineating and discovering biodiversity is critical for both taxonomic accuracy and conservation. In recent years, there has been a movement away from using exclusively morphological characters to delineate and describe taxa and an increase in the use of molecular markers to describe diversity or through integrative taxonomy, which employs traditional morphological characters, as well as genetic or other data. Tiger beetles are charismatic, of conservation concern, and much work has been done on the morphological delineation of species and subspecies, but few of these taxa have been tested with genetic analyses. In this study, we tested morphologically based taxonomic hypotheses of polymorphic tiger beetles in the Eunota circumpicta (LaFerté-Sénectère, 1841) species complex using multilocus genomic and mtDNA analyses. We find multiple cryptic species within the previous taxonomic concept of Eunota circumpicta, some of which were historically recognized as subspecies. We found that the mtDNA and genomic datasets did not identify the same taxonomic units and that the mtDNA was most at odds with all other genetic and morphological patterns. Overall, we describe new cryptic diversity, which raises important conservation concerns, and provide a working example for testing species and subspecies validity despite discordant data.

A near-complete genome assembly of Monochamus alternatus a major vector beetle of pinewood nematode.

The Japanese sawyer beetle, Monochamus alternatus, is not only one of the most important wood boring pest itself, but also a major vector of the invasive pinewood nematode (PWN), which is the causal agent of the devastative pine wilt disease (PWD) and threats the global pine forest. Here, we present a near-complete genome of M. alternatus at the chromosome level. The assembled genome was 792.05 Mb with contig N50 length of 55.99 Mb, which is the largest N50 size among the sequenced Coleoptera insects currently. 99.57% of sequence was anchored onto ten pseudochromosomes (one X-chromosome and nine autosomes), and the final genome harbored only 13 gaps. BUSCO evaluation revealed the presence of 99.0% of complete core genes. Thus, our genome assembly represented the highest-contiguity genome assembly as well as high completeness in insects so far. We identified 20,471 protein-coding genes, of which 20,070 (98.04%) were functionally annotated. The genome assembly of M. alternatus provides a valuable resource for exploring the evolution of the symbiosis between PWN and the vector insects.

The direction, timing and demography of Popillia japonica (Coleoptera) invasion reconstructed using complete mitochondrial genomes.

The Japanese beetle Popillia japonica is a pest insect that feeds on hundreds of species of wild and cultivated plants including important fruit, vegetable, and field crops. Native to Japan, the pest has invaded large areas of the USA, Canada, the Azores (Portugal), Italy, and Ticino (Switzerland), and it is considered a priority for control in the European Union. We determined the complete mitochondrial genome sequence in 86 individuals covering the entire distribution of the species. Phylogenetic analysis supports a major division between South Japan and Central/North Japan, with invasive samples coming from the latter. The origin of invasive USA samples is incompatible, in terms of the timing of the event, with a single introduction, with multiple Japanese lineages having been introduced and one accounting for most of the population expansion locally. The origin of the two invasive European populations is compatible with two different invasions followed by minimal differentiation locally. Population analyses provide the possibility to estimate the rate of sequence change from the data and to date major invasion events. Demographic analysis identifies a population expansion followed by a period of contraction prior to the invasion. The present study adds a time and demographic dimension to available reconstructions.